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Theory predicts the low temperature magnetic excitations in spin ices consist of deconfined mag- 
netic charges, or monopoles. A recent transverse-field (TF) muon spin rotation (/iSR) experiment [S 
T Bramwell et ai, Nature 461, 956 (2009)] reports results claiming to be consistent with the tem- 
perature and magnetic field dependence anticipated for monopole nucleation — the so-called second 
Wien effect. We demonstrate via a new series of /iSR experiments in Dy2Ti2 07 that such an effect 
is not observable in a TF jiSK experiment. Rather, as found in many highly frustrated magnetic 
materials, we observe spin fiuctuations which become temperature independent at low temperatures, 
behavior which dominates over any possible signature of thermally nucleated monopole excitations. 

PACS numbers: 75.50.Dd, 75.40. Cx, 75.40.Gb, 76.75.-fi 



Spin ices, such as Ho2Ti207 and Dy2Ti207, are top- 
ical highly frustrated magnetic systems which exhibit 
a gamut of very interesting phenomena [J, In 
(Ho/Dy)2Ti207, the Ho'^"'" and Dy'^^ magnetic ions re- 
side on the vertices of a pyrochlore lattice of corner- 
sharing tetrahedra. A large single ion anisotropy forces 
the moment to point strictly along local (111) crystalline 
axes, along the line which connects the centers of the two 
adjoining tetrahedra and their common vertex, making 
the moments classical "local" Ising spins. Since Ho''+ 
and Dy'^"*' carry a large magnetic moment of ^ 10 /iB, 
the dipolar interaction in these systems is ~ 1 K at 
nearest neighbor distance and of similar magnitude as 
the Curie- Weiss temperature ^cw 01 ■ The frustration 
in spin ices stems from the long-range nature of 

the magnetic dipolar interaction and its consequential 
"self-screening" (r is the distance between ions) 
As a result, spin ices are frustrated ferromagnets with 
low-energy states characterized by two spins "pointing 
in" and two spins "pointing out" on each tetrahedron - 
the 2-in/2-out rule which defines minimum energy spin 
configurations. These map onto the allowed proton con- 
figurations in water ice which obey the Bernal-Fowler 
ice-rules 0, hence the name spin ice. 

The dipolar spin ice model and its refinement [7,] 
yield an accurate microscopic quantitative description of 
the equilibrium thermodynamic properties of spin ices, 
both in zero and nonzero magnetic field. In contrast, the 
problem of the dynamical response of the moments in 
spin ices remains much less studied and understood. An 



exciting recent development in that direction is the real- 
ization that the "2-in/2-out" spin configurations may be 
described via a divergence-free coarse-grained magneti- 
zation density field [5, 8]. A thermal fluctuation causing 
the flip of an Ising spin from an "in" to an "out" direc- 
tion amounts to the creation of a nearest-neighbor pair 
of magnetization source and sink on the two adjoining 
tetrahedra or, in other words, to the nucleation of "mag- 
netic monopoles" out of the spin-ice-rule obeying vac- 
uum [5]. Particularly interesting is the observation that 
monopoles in dipolar spin ice interact via an emerging 
Coulomb potential which decays inversely proportional 
to the distance which separates them and are therefore 
deconfined [1] . A recent numerical study [§] provides evi- 
dence that the temperature dependence of the relaxation 
time determined in ac magnetic susceptibility measure- 
ments [lo| can be rationalized in terms of thermally ac- 
tivated monopoles, at least above 1 K. The wave vector 
dependence of the neutron scattering intensity suggests 
power law spin correlations, which are a prere quis ite for 
monopoles with effective Coulomb interactions Yet, 
perhaps the reported direct evidence for the presence of 
monopoles in spin ice and a determination of their effec- 
tive charge is the most intriguing recent result [l2| . 

In weak electrolytes, including water ice, characterized 
by a small dissociation rate constant K, the so-called 
second Wien effect describes the nonlinear increase of 
K under an applied electric field. In a recent paper fl3l , 
Bramwell and co-workers have drawn further on the anal- 
ogy between magnetic moments in spin ice and protons 
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in water ice [3. Using Onsager's accurate theory of the 
Wien effect |13j . Bramweh et al. put forward an elegant 
model to describe the dependence of the monopole nu- 
cleation rate, k(T, H) in spin ice on temperature T and 
external applied magnetic field H. They proposed that a 
measurement of k{T, H) could yield the monopole charge 



Q. To that effect, the authors of Ref. [l^ used ^SR in 
a transverse-field (TF) geometry to determine K(r, H) 
and extract a value Qexp ~ S/ie/ A, close to the value 
Qtheo ~ 4.6/iB/ A anticipated by theory |5|. 

In this Letter, we discuss how the weak TF /iSR ex- 
periment of Ref. [13], as a means to observe the second 
Wien effect in spin ice, was flawed in its conceptual de- 
sign and execution and incorrect in its theoretical inter- 
pretation of the muon spin depolarization rate. Monte 
Carlo calculations show that the internal magnetic field 
at the expected muon locations in Dy2Ti207 spin ice ma- 
terial is "large" 0.3 T) and has a broad distribution, 
preventing the observation of TF muon precession. We 
present evidence that the coherent muon precession in 
weak {H ^ 1 mT) TF seen in Ref. 12| originated rather 
from the sample holder and other parts of the sample en- 
vironment. In contrast, our zero-field /iSR results exhibit 
low-temperature muon spin relaxation which is temper- 
ature independent from 4 K down to 20 mK. 

Positive muons provide a pointlike real space magnetic 
probe averaging over the Brillouin zone, in contrast with 
magnetization measurements, which measure only the 
Q=0 response. In a ^SR experiment, essentially 100% 
spin-polarized positive muons are implanted in a mate- 
rial and precess in the local magnetic field B(r). The 
muons subsequently decay (with lifetime = 2.2 ^s) 
into a positron (emitted preferentially in the direction 
of the muon spin at the time of decay) and two unde- 
tected neutrinos. An asymmetry signal, obtained from 
the decay histograms of opposing positron detectors, rep- 
resents the projection of the muon spin polarization func- 
tion onto the axis defined by the detectors. Since the 
muons are created fully spin polarized (via the parity- 
violating weak decay of their parent pions), /iSR experi- 
ments may be performed in zero (ZF), longitudinal (LF) 
or transverse (TF) magnetic field. As in nuclear magnetic 
resonance (NMR), the depolarization results from both 
dynamic (Ti) and static {T2) processes. Further details 
describing /xSR are found elsewhere [14 1. 

TF-/xSR - In a TF experiment, an external magnetic 
field is applied perpendicular to the initial muon spin 
polarization direction. For arbitrary electronic spin fiuc- 
tuation rates, an approximate analytic form for the high 
TF-/iSR polarization function is given by a relaxation 
envelope multiplied by a cosine precession signal [3, [3] , 

P^,{t) ^ exp [-{IS^/v'^) (e-^^* - 1 + ut)] coa{ujt). (1) 

In Eq. (1), A = "(^,3 is the muon gyromagnetic ratio 
times the rms instantaneous internal magnetic field at 
the muon site B and v is the field fluctuation rate. In the 



fast fluctuation {v 3> A) regime, the envelope becomes 
exp(— A^i/z^) whereas, in the slow fluctuation (z/ <^ A) 
regime, the envelope reduces to exp(— A^f^/2), which is 
independent of v. Hence, the relaxation in a transverse 
field (TF) never takes an exponential form where the re- 
laxation rate l/Ti cx ly, as assumed by Bramwell et al. 
[l^ . As discussed further below, such behavior is only 
applicable to ZF and LF measurements. As the fluc- 
tuation rate of the electronic moments decreases and a 
significant spectral density develops near zero frequency, 
the local field at the muon site becomes the vector sum 
of the applied and internal fields. As we show, this net 
field in dipolar spin ice has a much larger rms value than 
the applied field H. Therefore the muon polarization 
function is not given by a cosine with a frequency cor- 
responding to the applied field but, instead, is rapidly 
damped to zero. 

From the discussion above, a crucial issue is whether 
the internal field distribution P{B{r)) in a spin ice has 
significant weight below the applied external field of 
H ^ 1 mT. As a first step to address this question, we 
use a loop algorithm 16| in Monte Carlo (MC) simula- 
tions of a reahstic microscopic model of Dy2Ti207 0| to 
calculate P{B{r)) at the most probable muon locations, 
as determined by density functional theory (DFT) cal- 
culations [l7| . We find, confirming naive expectations, 
that P{B{r)) is heavily populated for fields of several 
hundreds of millitesla. Such values are consistent with 
the estimate of ^ 0.5 T obtained by Lago et al. from 
LF- ^SR decoupling measurements in Dy2Ti207 [3]. In 
fact, for all four lowest energy potential muon stopping 
sites, we find vanishing P{B{r)) at \{B{r))\ 0. 

Figures [TJa) andlljb) show TF-^SR spectra measured 
in single crystals of Dy2Ti207 mounted using GE var- 
nish on an intrinsic GaAs plate (blue triangles), as well 
as results on a blank GaAs plate without the Dy2Ti207 
sample (red circles) with an external field of 2 mT in both 
cases. The crystals were grown using floating zone image 
furnace methods, with the speed of 4 mm/h. GaAs was 
chosen as it exhibits no precession at the muon Larmor 
frequency since all muons form muonium, a hydrogenlike 
muon-electron bound pair [l3|. As a result, the observed 
red precession signal in Figs. [D^a) and [TJ^b) is a purely 
instrumental background from muons which land else- 
where in the dilution refrigerator or in the silver (Ag) 
sample holder. The signal with the Dy2Ti207 sample on 
the GaAs plate (blue) is essentially identical to the back- 
ground contribution in Fig.[TKa) at long {t ^ 1 /is) times. 
This demonstrates that there is no long-lived muon pre- 
cession signal originating from the specimen, consistent 
with the aforementioned expectation based on DFT-MC 
calculations 17|. The small applied H = 2 mT is hence 
negligible compared to B{r). The small difference in 
Fig. [TJb) at early times (t ^ 1 /is) comes from the lon- 
gitudinal relaxation of the /iSR signal in the Dy2Ti207 
sample, visible here because the initial muon spin po- 
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larization is parallel to the positron detector axis. This 
difference is absent in Fig. [TJa) where, due to the chosen 
experimental geometry, the initial muon polarization is 
perpendicular to the detector axis. TF-/iSR spectra mea- 
sured with the Dy2Ti207 crystals mounted on a high pu- 
rity Ag plate for the same two experimental geometries 
are shown in black. Using a metal (e.g. silver) back- 
ing plate ensures the sample is in good thermal contact. 
Additionally, silver (taken alone) produces an essentially 
undamped TF-/iSR precession signal. The increased sig- 
nal amplitude relative to the signal using a GaAs back- 
ing reflects the contribution of muons landing in parts 
of the Ag backing not covered by the sample. Such TF- 
/iSR spectra measured using Ag are essentially identical 
to those reported in Ref. [12]. Note that in both in- 
stances, to account for history dependent effects, care 
was taken to cool the samples in zero field before ap- 
plying 2 mT at base temperature and taking data while 
warming. We note that the evolution of the TF depolar- 
ization rate tracks the magnetization (T^]. We speculate 
that a stray field, proportional to the dc magnetization, 
is generated within the Ag plate in the areas between the 
crystallites of Dy2Ti207. Muons landing in such a region 
undergo slow relaxation due to the inhomogeneous stray 
field, proportional to the applied field, consistent with 
observations of Bramwell et al. . 

The above experiments, performed in a dilution re- 
frigerator, require a cold finger sample holder to provide 
thermal contact to the sample. In contrast, by using a 
''He gas-flow cryostat operating at temperatures above 
T — 2 K, one can perform /xSR without complications 
due to a background signal by suspending the specimen 
on thin tape. The TF-^SR signal observed in such a 
"background free" apparatus is shown in Fig. [2][a) in 
H = 2 mT. The absence of a precession signal at T = 2 K 
confirms that any long lived precession signal seen at 
T ^ 2 K does not originate from Dy2Ti2 07. Rather, 
it is a background signal from the sample holder or cryo- 
stat. This constitutes the first of our two main results. 

ZF-/iSR - A more effective method for studying spin 
dynamics in magnetic systems is to measure the spin po- 
larization in the ZF/LF geometry At high tempera- 
tures, in the fast fluctuation regime, the spin lattice relax- 
ation rate is l/Ti = 2iS} jv in zero applied fleld As 
shown in Fig.|3l \jT\ increases as temperature decreases. 
This is due to the combined effect of changes in the size 
of A as the Dy^+ excited crystal electric field levels are 
depopulated and the slowing down of the Dy'^^ fluctu- 
ation rate v. The relaxation rate peaks at T « 50 K 
and drops below, entering the slow fluctuation regime, 
consistent with earlier studies on powder samples [l8| . 
Below T = 50 K the muon spin polarization exhibits a 
two component form, indicating that the local magnetic 
environment consists of a large quasistatic field [respon- 
sible for the rapid loss of polarization seen in Fig.[2l^b) at 




FIG. 1. (Color online) TF /iSR spectra measured at T = 
100 rtiK in Dy2Ti2 07 with two different counter geometries. 
In (a), the applied field H — 2 mT lies in the plane of the 
platelike coaligned mosaic of single crystals. In (b), H is 
perpendicular to the plates. In both cases H is along the 
crystallographic [100] direction. The positron detectors are 
indicated as back (B), forward (F), left (L) and right (R). 
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FIG. 2. (Color online) Background free (a) TF=2 mT and 
(b) ZF /iSR spectra measured in Dy2Ti2 07. Note the over- 
damped precession signal in (a). The positron detectors are 
indicated as back (B), forward (F), up (17) and down (D). 



2 K] coexisting with a fiuctuating field component (which 
gives the damping of the remaining polarization). The 
amplitude of the slowly relaxing component at T ^ 50 K 
is 1/3 of the initial polarization, as expected for a cubic 
material where 1 /3 of the muon polarization is on average 
parallel to B(r) Analogous decreases in Dy'^^ fluc- 
tuation rate over 4 orders of magnitude between 300 and 
8 K have been observed using zero field ''^Ti NQR 
and nuclear forward scattering [2l[. 

We observe a substantial 1/Ti ^ 1 /is^^ relaxation 
rate at T < 5 K, more than an order of magnitude 
above our detection limit (~ 10~^ MS"^)- This observa- 
tion contrasts dramatically with the activated behavior 
anticipated for magnetic monopoles Q and is startling 
given the highly Ising nature of the Dy^"*" spins, where 
large energy barriers against single spin flip processes 
separate quasidegenerate ice rules states [Ij. The ori- 
gin of this temperature independent relaxation is as yet 
unclear. We note that so-called persistent spin dynam- 
ics have been observed in a wide range of geometrically 
frustrated materials 2^, l23|. Oddly, the characteristic 
rare earth spin fluctuation rates extracted using vari- 
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FIG. 3. Muon spin relaxation rate in Dy2Ti207. Individ- 
ual /iSR spectra were analyzed using a phenomenological 
stretched exponential form Pz(t) x exp(— (t/Ti)'') commonly 
used to model glassy systems. The low temperature behavior 
of T]~^ is shown on an expanded linear scale in the inset. The 
exponent /3 drops monotonically from ~ 0.75 at 150 K to 0.4 
below 5 K. 



of unusual spin excitations in a model Ising system. It 
poses the challenge to comprehensively understand the 
microscopic mechanism(s) causing the temperature inde- 
pendent muon spin relaxation within the broader context 
of geometrically frustrated systems [2^ . 
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ous techniques differ dramatically in Dy2Ti207, as do 
the values when comparing the two isostructural com- 
pounds AaTiaOr (A=Dy,Ho) [13, The strong 
hyperfine interaction between the electronic and nuclear 
spin species should also not be neglected, particularly in 
the latter compound, as highlighted by the pronounced 
Schottky anomaly observed arising from nuclear contri- 
butions to the magnetic specific heat [2^. 

Dy-based compounds form a variety of model Ising 
systems (30j . However, many of them exhibit unex- 
pected dynamic spin fluctuations at low temperatures: 
(1) the single molecular magnetic system [DyPc2]°, char- 
acterized by a doubly degenerate ground state and large 
magnetic anisotropy, exhibits a tunnelling regime [3l|, 
suggesting such temperature independent behavior is a 
more pervasive phenomenon [s^ : (2) the geometrically 
frustrated Ising antiferromagnet Dysprosium Aluminium 
Garnet, where marked changes in the characteristic re- 
laxation times over several orders of magnitude have been 
reported as a function of applied fields ^30,]; (3) even the 
archetypal dilute Ising dipolar ferromagnet Dysprosium 
Ethyl Sulphate, which exhibits unex pec tedly high relax- 
ation rates within the ordered state 
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Despite a seemingly compelling argument for spin dy- 
namics caused by monopoles 0, additional spin re- 
laxation processes dominate the behavior observable us- 
ing ZF-/.tSR in Dy2Ti207 spin ice. In Heisenberg spin 
systems, geo metrical frustration may lead to spin liquid 
behavior [3J] where spin fluctuations persist to absolute 
zero. Understanding the low-temperature dynamics in 
the Ising Dy2Ti207 system will require the construction 
of an effective low-energy Hamiltonian containing non- 
Ising terms [s^. The present work clarifies the nature 
of static and dynamic contributions to the internal mag- 
netic field as probed using /iSR and describes evidence 
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SUPPLEMENTARY MATERIAL FOR: 
MAGNETIC EXCITATIONS WITHOUT 
MONOPOLE SIGNATURES USING MUON SPIN 
ROTATION 

This document is a supplement to our main article, in 
which we provide some of the details of the density func- 
tional theory (DFT) and Monte Carlo (MC) calculations 
which were carried out to determine the most probable 
muon location sites in the Dy2Ti207 spin ice compound, 
as well as the associated internal magnetic field distribu- 
tion. 



DFT calculations 

To determine the quasistatic dipolar magnetic field 
which a positive muon experiences in the insulating 
Dy2Ti2 07 spin ice compound, one must first determine 
the most probable local minimum energy locations a 
muon occupies after coming to rest upon high (4.1 MeV 
kinetic) energy implantation in the sample. The stan- 
dard expectation is that a positive muon (a light version 
of proton) forms a (hydrogenlike) bond with an oxygen 
at a typical distance of 1 A. To address this question 
more quantitively, we determined those low-energy muon 
locations (ground state and a few metastable states) us- 
ing density functional theory (DFT) calculations using 
the VASP package 0. In those calculations, one single 
positive muon was placed inside an 88- atom Dy2Ti2 07 
supercell with periodic boundary conditions, which was 
neutralized by applying a uniform negative background 
charge. The valence wavefunctions were expanded in 
a plane-wave basis with a cutoff energy of 400 eV. A 
2x2x2 grid was used for the fc-point sampling of 
the Brillouin zone. Perdew-Burke-Ernzerhof exchange- 
correlation functionals [l[ were used and all atoms in 
the cell were relaxed to minimize the Feynman-Hellmann 
forces to below 0.01 ev/A. The ground state as well as 
three other metastable higher energy states were identi- 
fied, with their energy measured from the ground state 
labeled in Fig. The determined ground state as well 
as the identified metastable higher energy states of muon 



locations were taken as the most likely candidate for the 
muon locations in Dy2Ti207. 
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FIG. 4. (Color online). Probability histograms of the internal 
magnetic fields at candidate muon sites are sampled by Monte 
Carlo (MC) simulations at 200 mK. The insert shows the 
time-dependence of the muon spin polarization for both the 
muon location ground state (curve labeled ground) and the 
average of the four lowest energy states. 



MC calculations 

The internal magnetic fields at the muon sites due to 
the static Dy^+ magnetic moments were next calculated. 
To determine this field, we performed Monte Carlo (MC) 
simulations to sample the magnetic moment configura- 
tions at various temperatures. For these calculations, any 
local structural distortions of the original crystal caused 
by the presence of the muon were ignored in the MC. 
We considered the same extended dipolar spin ice model 
of Dy2Ti207 as in Ref. [1] with exchange interactions 
up to third nearest neighbors and long-range magneto- 
static dipole-dipole interactions. We employed a system 
size of 432 spins with periodic boundary conditions. The 
dipole-dipole interactions among the Dy'^'^ magnetic mo- 
ments as well as the local internal magnetic field at the 
muon sites determined by DFT were calculated using the 
Ewald summation technique. A loop update algorithm 
0] was employed to overcome the energy barriers of the 
single-spin-fiip updates at temperatures below the Schot- 
tky peak in the specific heat of Dy2Ti207, which occurs 
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at a temperature T^ak ~ 1-1 K. At least 1000 MC steps 
with loop updates [4| per spin were used in equilibrating 
the system. The sampling of the local dipolar magnetic 
field was done using 360,000 spin configurations. We as- 
sume that the time scale for the spin dynamics is, deep 
in the spin ice regime at low temperature (T <C Tpeak), 
much longer than the muon lifetime and that the muons 
implanted homogeneously in the sample, one by one, ex- 
perience static fields from these spin configurations with 
equal probability. No muon-Dy'^"'" transferred hyperfine 
field was taken into account in the present calculations. 
The probability histograms of the internal magnetic field 
at 200 mK are plotted in the main panel of Fig. 21 The 
multiple peaked structures are a combined consequence 
of (i) of the discreteness of the Ising spins and (ii) the 
fact that spin ices are not completely random as in an 
Edwards- Anderson spin glass, but rather, obey the two- 
in/two-out ice rules with the consequential dipolar cor- 
relations 0. 

One observes that the probability density for magnetic 
fields of order of hundreds of millitesla or higher is large 
at the lowest energy candidate muon sites determined by 
DFT. As an independent check, we determined the in- 
ternal magnetic field at the surface of a spherical shell 
of radius of 1 A centered on oxygen locations and con- 
firmed that such a magnitude of internal magnetic field 
is indeed typical. We also note that there is vanishing 
probability of B{r) < 10 mT for all muon locations. The 
inset of Fig. 2] shows the zero field (ZF) muon spin depo- 
larization associated with the muon ground state (curve 
labeled "ground") as well as the average relaxation for 
the ground state and the three lowest metastable states 
(curve labeled "4 site average") with initial muon po- 
larization along [100] (see main text). The fast oscilla- 
tions of the muon polarization for short time (t < 0.5 /is) 
with its rapid decay and asymptotic approach to a 1/3 
value is again a reflection of the typical magnitude of the 
magnetic field a muon experiences, which is of order 100 
mT. Interestingly, the experimental /iSR data do not find 
such "front end" oscillations, which should be observable 
given the time binning used in the data acquisition and 
the instrumental (electronic) timing resolution. This is 
possibly a result of the slight randomization of the local 
crystal field Hamiltonian acting at the Dy'^"'" ions neigh- 
boring the implanted muon, and hence the corresponding 



weakly randomized magnitude of the Dy^"*" magnetic mo- 
ments Q, as well as their direction and which we have 
ignored. Another randomization which we have ignored 
is that of the transferred hyperfine field. Finally, and 
most notably, these DFT-MC calculations neglect the ef- 
fect of Dy'^^ spin fluctuations on the timescale of the 
muon lifetime — perhaps the most signiflcant new result 
of our experimental study and reported in the main part 
of the paper. Indeed, as observed experimentally, these 
have perhaps the most dramatic effect on the muon spin 
depolarization function. 

While all these effects which go beyond the approxima- 
tions used here should ultimately be taken into account 
in future calculations, the results of such improved calcu- 
lations will not affect our main conclusion. Namely, the 
conflguration of Dy'^^ magnetic moments in Dy2Ti207 
are sufficiently random, despite the ice-rule, that the like- 
lihood that a finite percentage of muons land at sites with 
a local magnetic field negligible compared to the applied 
external field of 2 mT is for all practical purposes zero. 
Ultimately, the experimental results reported in the main 
part of the paper corroborate this assertion: no weak TF 
precession signal whose origin is the Dy2Ti207 sample 
can be positively identified, confirming the above DFT- 
MC calculations that there are no density of (muon oc- 
cupation) states with local field of the order or less than 
the applied magnetic fields of order of 2 mT. 
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